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C\| ■ Abstract 

A study of energy behavior of the pion spectra and interferometry scales 
is carried out for the top SPS, RHIC and for LHC energies within the hy- 
drokinetic approach. The main mechanisms that lead to the paradoxical, at 
^ ' first sight, dependence of the interferometry scales with an energy growth, in 

particular, a decrease R ou t/Rside ratio, are exposed. The hydrokinetic pre- 
dictions for the HBT radii at LHC energies are compared with the recent 
results of the ALICE experiment. 
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1 Introduction 



In anticipation of start of the LHC experiments there were presented many 
■ different theoretical views on the underlying physics at so large collision 

energies as well as numerous predictions for observables - a big collection 
of the predictions for heavy ion collision was assembled in [T]. The similar 
took place on the eve of the RHIC experiments. At that "pre-RHIC" time 
the future results on the correlation femtoscopy of particles, that is the topic 
of this entry, were expecting with great interest. One of the reason was a 
hope to find the interferometry signature of the quark-gluon plasma (QGP). 
^ ! A very large value of the ratio of the two transverse interferometry radii, 

Rout to R S ide, was predicted as a signal of the QGP formation While 
the R S ide radius is associated with the transverse homogeneity length [3], the 
R ou t includes besides that also additional contributions, in particular, the 
one which is related to a duration of the pion emission. Since the lifetime of 
the systems obviously should grow with collision energy, if it is accompanied 
by an increase of the initial energy density and/or by a softening of the 



1 The entry is based on the talks given at the Sixth Workshop on Particle Cor- 
relations and Femtoscopy, BITP, Kiev, September 14 - 18, 2010 and GSI/EMMI 
Seminar, January 14, 2011 
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equation of state due to phase transition between hadron matter and quark- 
gluon plasma (QGP), the duration of pion emission should also grow with 
energy and so R out /R S id e ratio could increase. 

The RHIC experiments brought an unexpected result: the ratio 
Rout/Rside ~ 1 is similar or even smaller than at SPS. The another surprise 
was the absolute values of the radii. Naively it was expected that when the 
energy of colliding nuclei increases, the pion interferometry volume Vi nt - 
product of the interferometry radii in three orthogonal directions - will rise 
at the same maximal centrality for Pb+Pb and Au+Au collisions just pro- 
portionally to . However, when experiments at RHIC starts, an increase 
of the interferometry volume with energy turn out to be essentially smaller 
then if the proportionality law takes place. Both these unexpected results 
were called the RHIC HBT puzzle @]. 

During a long period this puzzle was not solved in hydrodynamic/hybrid 
models of A+A collisions which reproduce good the single particle trans- 
verse spectra and its axial anisotropy in non-central collisions described by 
the v 2 coefficients. Only a few years ago the main factors which allow one 
to describe simultaneously the spectra and femtoscopic scales at RHIC be- 
come clear. They are [5]- [10]: a relatively hard equation of state because of 
crossover transition (instead of the 1st order one) between quark-gluon and 
hadron phases and due to nonequilibrium composition of hadronic matter, 
presence of prethermal transverse flows and their anisotropy developed to 
thermalization time, an 'additional portion' of the transverse flows owing to 
the shear viscosity effect and fluctuation of initial conditions. An account of 
these factors gives the possibility to describe good the pion and kaon spectra 
together with the femtoscopy data of RHIC within realistic freeze-out picture 
with a gradual decay of nonequilibrium fluid into observed particles |11| . 

Now, when the heavy ion experiments at LHC already starts, and the 
ALICE Collaboration published the first results on the femtoscopy in A+A 
collisions at y/s = 2.76 TeV [12], the main question is whether an understand- 
ing of the physics responsible for the space-time matter evolution in Au+Au 
collisions at RHIC can be extrapolated to the LHC energies, or some new 
"LHC HBT puzzle" is already apprehended just as it happened in the way 
from SPS to RHIC energies. In this note we describe the physical mechanisms 
responsible for the peculiarities of energy dependence of the interferometry 
radii and therefore solving the RHIC HBT puzzle, present the quantitative 
predictions given for LHC within hydrokinetic model earlier [13] . compare 
them with the recent ALICE LHC results and make the corresponding infer- 
ence. 
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2 Hydro-kinetic approach to A+A collisions 



Let us briefly describe the main features of the HKM [THE]. It incorporates 
hydrodynamical expansion of the systems formed in A+A collisions and their 
dynamical decoupling described by escape probabilities. 

Initial conditions — Our results are all related to the central rapidity slice 
where we use the boost-invariant Bjorken-like initial condition. We consider 
the proper time of thermalization of quark-gluon matter to be To = 1 fm/c, 
at present there is no theoretical arguments permitting smaller value. The 
initial energy density in the transverse plane is supposed to be Glauber-like 
|15j . i.e. is proportional to the participant nucleon density for Pb+Pb (SPS) 
and Au+Au (RHIC, LHC) collisions with zero impact parameter. The height 
of the distribution - the maximal initial energy density - e(r = 0) = eo is the 
fitting parameter. From analysis of pion transverse spectra we choose it for 
the top SPS energy to be e = 9 GeV/fm 3 ((e) = 6.4 GeV/fm 3 ), for the 
top RHIC energy e = 16.5 GeV/fm 3 ((e) = 11.6 GeV/fm 3 ). The brackets 
< ... > correspond to mean value over the distribution associated with the 
Glauber transverse profile. We also demonstrate results at eo = 40 GeV/fm 3 
and eo = 60 GeV/fm 3 . In hydrokinetic model eo = 40 GeV/fm 3 correspond to 
multiplicity of charged particles dN c h/drj 1500. We suppose that soon after 
thermalization the matter created in A+A collision at energies considered is 
in the quark gluon plasma (QGP) state. 

At the time of thermalization, tq = 1 fm/c, the system already has de- 
veloped collective transverse velocities 0E]. The initial transverse rapidity 
profile is supposed to be linear in radius Tt- 

y T = (1) 



where a is the second fitting parameter and Rt = a/< r\ >. Note that 
the fitting parameter a should absorbs also a positive correction for under- 
estimated resulting transverse flow since in this work we did not account in 
direct way for the viscosity effects [16] neither at QGP stage nor at hadronic 
one. In formalism of HKM [6] the viscosity effects at hadronic stage are 
incorporated in the mechanisms of the back reaction of particle emission on 
hydrodynamic evolution which we ignore in current calculations. Since the 
corrections to transverse flows which depend on unknown viscosity coeffi- 
cients are unknown, we use fitting parameter a to describe the "additional 
unknown portion" of flows, caused both factors: by a developing of the pre- 
thermal flows and the viscosity effects in quark-gluon plasma. The best fits 
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of the pion transverse spectra at SPS and RHIC are provided at a = 0.194 
((v T ) = 0.178) for SPS energies and a = 0.28 ((v T ) = 0.25) for RHIC ones. 
The latter value we use also for LHC energies aiming to analyze just influence 
of energy density increase. 

Equation of state — Following to Ref. [6] we use at high temperatures 
the EoS [T7] adjusted to the QCD lattice data with the baryonic chemical 
potential Hb = and matched with chemically equilibrated multi-component 
hadron resonance gas at T = 175 MeV. Such an EoS could be a good approx- 
imation for the RHIC and LHC energies; as for the SPS energies we utilize 
it just to demonstrate the energy dependent mechanism of formation of the 
space-time scales 0. We suppose the chemical freeze-out for the hadron gas 
at T c h = 165 MeV [19J. It guarantees us the correct particle number ratios 
for all quasi-stable particles (here we calculate only pion observables) at least 
for RHIC. Below T c h a composition of the hadron gas is changed only due to 
resonance decays into expanding fluid. We include 359 hadron states made 
of u, d, s quarks with masses up to 2.6 GeV. The EoS in this non chemically 
equilibrated system depends now on particle number densities Hi of all the 
359 particle species i: p = p(e, {fii}). Since the energy densities in expand- 
ing system do not directly correlate with resonance decays, all the variables 
in the EoS depend on space-time points and so an evaluation of the EoS is 
incorporated in the hydrodynamic code. We calculate the EoS below T c h in 
the Boltzmann approximation of ideal multi-component hadron gas. 

Evolution — At the temperatures higher than T c ^ the hydrodynamic evo- 
lution is related to the quark-gluon and hadron phases which are in chemical 
equilibrium with zero baryonic chemical potential. The evolution is described 
by the conservation law for the energy-momentum tensor of perfect fluid: 

d v T^[x) = (2) 

At T < T c h=lQ5 MeV the system evolves as non chemically equilibrated 
hadronic gas. The concept of the chemical freeze-out implies that afterwards 
only elastic collisions and resonance decays take place because of relatively 
small densities allied with a fast rate of expansion at the last stage. Thus, in 
addition to §2§ , the equations accounting for the particle number conservation 
and resonance decays are added. If one neglects the thermal motion of heavy 
resonances the equations for particle densities rii(x) take the form: 

d^(n i (x)u tl (x)) = -Yim{x) + ^ /;,,!' ,//,(./•) (3) 

j 

2 a good description of the spectra and HBT radii at the SPS energies with realistic 
EoS within hydrokinetic model is presented in Ref. 1 18) 
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where bij = BijNij denote the average number of i-th particles coming from 
arbitrary decay of j-th resonance, B^ = Tij/Tj ttot is branching ratio, is 
a number of i-th particles produced in j — > i decay channel. We also can 
account for recombination in the processes of resonance decays into expand- 
ing medium just by utilizing the effective decay width T^e// = 7IV We use 
7 = 0.75 supposing thus that near 30% of resonances are recombining during 
the evolution. All the equations and 359 equations ([3]) are solving simul- 
taneously with calculation of the EoS, p(x) = p(e(x), {rii(x)}), at each point 
x. 

System's decoupling and spectra formation — During the matter evolu- 
tion, in fact, at T < T c ^, hadrons continuously leave the system. Such a 
process is described by means of the emission function S(x,p) which is ex- 
pressed for pions through the gain term, G n (x,p), in Boltzmann equations 

00 

and the escape probabilities V n (x,p) = exp(— j dsR n+ h{s, r + ^(s — t),p)): 

S v (x,p) = G n (x,p)V n (x,p) [THE]. For pion emission in relaxation time ap- 
proximation G n ~ f n Rn+h + Gh->tt where f n (x,p) is the pion Bose-Einstein 
phase-space distribution, R n+ h(x,p) is the total collision rate of the pion, 
carrying momentum p, with all the hadrons h in the system in a vicinity of 
point x, the term Gh^-k describes an inflow of the pions into phase-space 
point (x,p) due to the resonance decays. It is calculated according to the 
kinematics of decays with simplification that the spectral function of the res- 
onance H is 5(p 2 — (mH) 2 )- The cross-sections in the hadronic gas, that deter- 
mine via the collision rate Rn+h the escape probabilities V(x,p) and emission 
function S(x,p), are calculated in accordance with the UrQMD method [20 J . 
The spectra and correlation functions are found from the emission function 
S in the standard way (see, e.g., |14|). 

3 Results and conclusions 

The pion emission function per unit (central) rapidity, integrated over az- 
imuthal angular and transverse momenta, is presented in Fig. 1 for the top 
SPS, RHIC and LHC energies as a function of transverse radius r and proper 
time r. The two fitting parameters and (vt) are fixed as discussed above, eo 
is also marked in figures. The pion transverse momentum spectrum, its slope 
as well as the absolute value, and the interferometry radii, including R out to 
R S ide ratio, are in a good agreement with the experimental data both for the 
top SPS and RHIC energies. 
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Figure 1: The p^-integrated emission functions of negative pions for the 
top SPS, RHIC and LHC energies (top); the interferometry radii (middle) 
R ut/Rside ratio and transverse momentum spectra (bottom) of negative pi- 
ons at different energy densities, all calculated in HKM model. The experi- 
mental data are taken from CERES [27] and NA-49 Collaborations [281 ES] 
(SPS CERN), STAR [301 E] and PHENIX [321 [33] Collaborations (RHIC 
BNL) and ALICE Collaboration (LHC, CERN) [12]. 

As one can see particle emission lasts a total lifetime of the fireballs; in the 
cental part, r « 0, the duration is half of the lifetime. Nevertheless, accord- 
ing to the results [HI [21], the Landau/Cooper- Frye presentation of sudden 
freeze-out could be applied in a generalized form accounting for momentum 
dependence of the freeze-out hypersurface <J p (x); now cr p (x) corresponds to 
the maximum of emission function S(t a (r,p),r,p) at fixed momentum p in 
an appropriate region of r. This finding allows one to keep in mind the known 
results based on the Cooper-Frye formalism, applying them to a surface of 
the maximal emission for given p. Then the typical features of the energy 
dependence can be understood as follows. The inverse of the spectra slopes, 
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T e ff, grows with energy, since as one sees from the emission functions, the 
duration of expansion increases with initial energy density and, therefore, the 
pressure gradient driven fluid elements gets more transverse collective veloc- 
ities vt when reach a decoupling energy densities. Therefore the blue shift 
of the spectra becomes stronger. A rise of the transverse collective flow with 
energy leads to some compensation of an increase of R S id e '- qualitatively the 
homogeneity length at decoupling stage is R S id e = Rceom/ a/1 + (v^)m T /2T, 
(see, e.g., |22|). So, despite an significant increase of the transverse emission 
region, Rceom, seen in Fig.l, a magnification of collective flow partially com- 
pensates this. It leads to only a moderate increase of the R S id e with energy. 
Since the temperatures in the regions of the maximal emission decrease very 
slowly when initial energy density grows (e.g., the temperatures for SPS, 
RHIC and LHC are cor respond ingly 0.105, 0.103 and 0.95 MeV for p T = 0.3 
GeV/c ) the Ri on g ~ r^/T/my |23j grows proportionally to an increase of the 
proper time associated with the hypersurface a PT (x) of maximal emission. As 
we see from Fig. 1 this time grows quite moderate with the collision energy. 

A non trivial result concerns the energy behavior of the R ou t/Rside ra- 
tio. It slowly drops when energy grows and apparently is saturated at fairly 
high energies at the value close to unity (Fig.l). To clarify the physical rea- 
son of it let us make a simple half-quantitative analysis. As one can see 
in Fig. 1, the hypersurface of the maximal emission can be approximated 
as consisting of two parts: the "volume" emission (V) at r « const and 
"surface" emission (S). A similar picture within the Cooper- Frye prescrip- 
tion, which generalizes the blast-wave model |24] by means of including of 
the surface emission has been considered in Ref. |25| . If the hypersurface 
of maximal emission f(r) is double-valued function, as in our case, then at 
some transverse momentum p? the transverse spectra and HBT radii will 
be formed mostly by the two contributions from the different regions with 
the homogeneity lengths A^y = a/ < (Ar^) 2 > (i = side, out) at the V- 
hypersurface and with the homogeneity lengths Aj )( g at the S-hypersurface. 
Similar to Ref. [22], one can apply at m^/T ^> 1 the saddle point method 
when calculate the single and two particle spectra using the boost-invariant 
measures fxy = da^p^ = r(r)rdrd(j)dr](m T cosh(?7 — y) — p T ^^- cos((p — a)) 

and = da^p^ = r(T)rdrd(j)dri(—mTCOsh.(ri — y)^jp- +PtCOs(</> — a)) for 
V- and S- parts of freeze-out hypersurface correspondingly (here r\ and y are 
space-time and particle pair rapidities, the similar correspondence is for an- 
gles and a, also note that ^ > [HI EI]). Then one can write, ignoring 
for simplicity the interference (cross-terms) between the surface and volume 
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contributions, 



H-side — c V^side,V + c S^side,S (^) 
o 2 \ 2 2 \ 2 / df , 2 / \ 

where the coefficients Cy + c| < 1 and we take into account that at p°/T 1 
for pions /3 01lt = p ou t/p° ~ 1. All homogeneity lengths depends on mean 
transverse momentum of the pion pairs p?. The slope j- in the region of 
homogeneity expresses the strength of r — r correlations between the space 
and time points of particle emission at the S-hypersurface r(r). The picture 
of emission in Fig. 1 shows that when the energy grows the correlations 
between the time and radial points of the emission becomes positive, ^ > 0, 
and they increase with energy density. The positivity is caused by the initial 
radial flows [5] tt r (ro), which are developed at the pre-thermal stage, and 
the strengthening of the r — r correlations happens because the non-central 
zth fluid elements, which produce after their expansion the surface emission, 
need more time Ti(eo) to reach the decoupling density if they initially have 
higher energy density 6q- (Let us characterized this effect by the parameter 
0). Then the fluid elements before their decays run up to larger 
radial freeze-out position rf if a is the average Lorentz-invariant acceleration 
of those fluid elements during the system expansion, then roughly for ith fluid 
elements which decays at time n we have at arj ^> 1: r^Tj) ~ rj(ro) + Tj + 
(m£(to) — l)/a. Then the level of r — r correlations within the homogeneous 
freeze-out "surface" region, which is formed by the expanding matter that 
initially at r occupies the region between the transversal radii T\{tq) and 
r 2 (r ) > ri(r ), is 

dr_ _ njn) - r 2 (r 2 ) _ ^ _ R_ 
dr T\ — t 2 e K 

and, therefore, the strength of r — r correlations grows with energy: ^ -> 1. 
Note that here we account for r 2 — T\ ~ K(eo(r 2 (ro)) — eo(ri(ro))) and that 
de °^ ~ — ^ where eo = eo( r = 0) and R is radius of nuclear. As a result 
the second S-term in Eq. (jSj) tends to zero at large eo , reducing, therefore, 
the R ou t/Rside ratio. In particular, if \ 2 sideV 3> A^ i(Je>s then, accounting for a 
similarity of the volume emission in our approximation and in the blast wave 
model, where as known \ S ide,v ~ ^outy, one can get: jj***- ~ 1+const--^ — > 1 
at eo — > oo. It is worthy of note that also measure fis tends to zero when 
^f- — > 1 that again reduces the surface contribution to side— and out— radii 
at large p T . 
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The presented qualitative analysis demonstrates the main mechanisms re- 
sponsible for the non-trivial behavior of R out to R S ide ratio exposed in HKM 
calculations, see Fig.l (bottom). The very recent first LHC data for Pb+Pb 
collisions presented by the ALICE Collaboration [12] conform, in fact, the 
discussed above physical picture of space-time evolution responsible for for- 
mation of the HBT radii and R out to R S id e ratio, see Fig.l. The transverse 
femtoscopy scales, predicted for the charged multiplicity dN c h/dr)=1500 in 
HKM at the initial energy density €q = 40 GeV/fm 3 , are quite close to the 
experimental data associated with dN c h/dr] ~ 1600 at the collision energy 
y/s = 2.76 TeV. As for the longitudinal HBT radius, Ri ong , it is underes- 
timated in HKM by around 20%. As the result, HKM gives smaller inter- 
ferometry volume than is observed at LHC. The reason could be that HKM 
describes a gradual decay of the system which evolves hydrodynamically until 
fairly large times. It is known [26] that at the isentropic and chemically frozen 
hydrodynamic evolution the interferometry volume increases quite moderate 
with initial energy density growth in collisions of the same/similar nucleus. 
The RHIC results support such a theoretical view (see solid line in Fig. 2), 
while the ALICE Collaboration observes a significant increase of the interfer- 
ometry volume at LHC. One should change, thus, the global fit of Vi n t(dN/drj) 
for A+A collisions for steeper slope (upper dash line). However, no one lin- 
ear fit cannot be extrapolated to Vi n t(dN/dr}) -dependence discovered by the 
ALICE Collaboration in p+p collisions [31] (bottom dashed line in Fig. 2). 
Could one call these two observed peculiarities as the '"LHC HBT puzzle 1 "? 
On our opinion, at least qualitavely, it is not puzzling. An essential growth 
of the interferometry volume in Pb+Pb collisons at the first LHC energy can 
be conditioned by an increase of the duration of the last very non-equilibrium 
stage of the matter evolution which cannot be considered on the hydrody- 
namic basis and one should use hadronic cascade models like UrQMD. At 
such late stage the results obtained in |26j for isentropic and chemically 
frousen evolution are violated. As for the different linear Vi n t(dN / drj) depen- 
dence in A+A and p+p collisions, the interferometry volume depends not 
only on multiplicity but also on intial size of colliding systems [26J. There- 
fore, qualitavely, we see no puzzle in the newest HBT results obtained at 
LHC in Pb+Pb and p+p collisions, but the final concusion can be done only 
after detailed quantitative analysis. 

Summary — We conclude that energy behavior of the pion interferometry 
scales can be understood at the same hydrokinetic basis as for the SPS and 
RHIC energies supplemented by hadronic cascade model at the latest stage 
of the evolution. In this approach the EoS accounts for a crossover transition 
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between quark-gluon and hadron matters at high collision energies and non- 
equilibrated expansion of the hadron-resonance gas at the later stage. 
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Figure 2: Illustration of multiplicity dependence of the pion inteferometry 
volume on charged particle multiplicity for central heavy ion collisions at 
AGS, SPS, RHIC and LHC energies and comparison with the results in p+p 
collisions (bottom, left). All the HBT radii are taken at the pion transverse 
momentum = 0.3 GeV. For A+A collision the data are taken from Fig. 4 
of Ref. [12j (see all details there), for p+p collisions the points for pt = 0.3 
GeV are interpolated from the results of Ref. [3l]. The solid line corresponds 
to linear fit of \^ nt ((iA r c / l /o?r7)-dependence only for the top SPS and the RHIC 
energies; upper dashed line is fit for all the A+A energies including the newest 
LHC point at s/s = 2.76 TeV; the bottom dashed line is the linear fit for 
ALICE LHC results for p+p collisions with energies 0.9 and 2.76 TeV. 



The HKM alows one to treat correctly the process of particle emission 
from expanding fireball, that is not sudden and lasts about system's lifetime. 
Also it takes into account the prethermal formation of transverse flows. Then 
the main mechanisms that lead to the paradoxical behavior of the interfer- 
ometry scales find a natural explanation. In particular, a slow decrease and 
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apparent saturation of R ou t/Rside ratio around unity at high energy happens 
due to a strengthening of positive correlations between space and time posi- 
tions of pions emitted at the radial periphery of the system. Such an effect 
is a consequence of the two factors accompanying an increase of collision 
energy: a developing of the pre-thermal collective transverse flows and an 
increase of initial energy density in the fireball. The prediction of the HKM 
for LHC energies are quite close to the first experimental data in Pb+Pb 
collisions at LHC. 
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